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We present a functional renormalization group analysis of a quantum critical point in a two-dimensional 
metal involving Fermi surface reconstruction due to the onset of spin density wave order. The critical theory is 
controlled by a fixed point in which the order parameter and fermionic quasiparticles are strongly coupled, and 
acquire spectral functions with a common dynamic critical exponent. We obtain results for critical exponents, 
and for the variation in the quasiparticle spectral weight around the Fermi surface. 
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Introduction. Quantum phase transitions between two 
Fermi liquids, one of which spontaneously breaks transla- 
tional symmetry and so reconstructs its Fermi surface, have 
been of long standing theoretical and experimental interest. 
Important new examples of experimental realizations have 
emerged in the past few years IBlSl^ and so a full theoretical 
understanding is of some urgency. However, despite several 
decades of theoretical work, key questions remain open espe- 
cially in the important case of two spatial dimensions. 

Early theories |4-^ for such quantum phase transitions fo- 
cused on effective models for the quantum fluctuations of the 
order parameter, while treating the Fermi surface reconstruc- 
tion as an ancillary phenomenon. However, it has since be- 
come clear ifTOl that such an approach is inadequate, and the 
Fermi surface excitations are primary actors in the critical 
theory. Ref. 11 postulated a critical theory for Fermi sur- 
face reconstruction, in which the Fermi surface excitations 
and the bosonic order parameter were equally important and 
both acquired anomalous dimensions. These excitations were 
strongly coupled to each other by a 'Yukawa' coupling of uni- 
versal strength, and their correlators scaled with a common 
dynamic critical exponent, z. Explicit computations were per- 
formed in the context of a 1 /A/^ expansion (where N is the num- 
ber of fermion flavors) for the onset of spin density wave or- 
der, but no such critical theory appeared at the two-loop level. 
Indeed, it was pointed out that at higher loops iFTTHTSll there is 
a breakdown of the 1 /N expansion, and so it remained unclear 
whether the postulated fixed point existed. 

Here we will address the problem of Fermi surface recon- 
struction at the onset of spin density wave order by an analysis 
based on a formally exact functional renormalization group 
approach In an approximation where we truncate the 

flow equations to a set of discrete points on the Fermi surface 
and keep the shape of the Fermi surface fixed, we establish the 
existence of a fixed point with the scaling structure postulated 
in Ref. [TTl describing the quantum phase transition between 
two Fermi liquids. A significant feature of this method is that 
it ties the parameters controlling the order parameter fluctua- 
tions to those associated with the fermion excitations, and this 
is important for a proper description of the scaling structure. 
We present numerical results for the critical exponents of the 
boson and fermion spectral functions, and for the variation in 
the fermionic quasiparticle residue around the Fermi surface. 



Model. Our computation will be carried in the context of 
the 'spin-fermion' model of antiferromagnetic fluctuations in 
a Fermi liquid |9 |. This involves a spin density wave order 
parameter at wavevector K = (tt, tt) coupled to fermions W 
moving on a square lattice. The analytic analyses have fo- 
cused on the vicinity of the 'hot spots' on the Fermi surface: 
these are the 8 points on the Fermi surface which can gener- 
ically be connected to each other by K. The fermion dis- 
persions were linearized and truncated around the hot spots. 
However a complete analysis requires that we avoid the spu- 
rious singularities associated with truncated Fermi surfaces, 
and deal only with continuous Fermi surfaces. Here, we will 
choose the Fermi surface configurations of a recent analysis 
|[T6ll which allowed Monte Carlo studies without a sign prob- 
lem. This paves the way for an eventual comparison of our 
renormalization group results with Monte Carlo. Our present 
method applies also to general Fermi surfaces, and provides 
access to real-time spectral functions which are not easily ob- 
tainable from imaginary-time Monte Carlo. 

The model of Ref. 1 16] contains fermions in two bands, ^q,, 
a = 1,2 (although our present method can also be applied to 
single band models) coupled to in the eflTective action 

+ j^^k-q) [q^ + r) $(q) (1) 

Jk,q 

where represents integrals over spatial momenta k = (kx,ky) 
over the Brillouin zone, and over frequencies ^o- The fermion 
spinors are defined by ^a(k) = {^a^ik) 4^a,i(k)), o; = 1,2. 
We already introduce here the cutoflT A along which we later 
integrate our renormalization group flow toward A ^ 0. With 
A = Auv we have the bare lattice action. The boson quadratic 
terms consists of the control parameter r and a spatial gradient 
squared to account for spatial variations of the order parameter 
field 0. The quantum dynamics of will be generated in the 
RG flow; putting a term into Eq. ([TJ) does not change our 
results. The fermion dispersions for nearest-neighbor hopping 
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FIG. 1: (Color online) Reconstructing Fermi surfaces (^k,i = 0, 
black-dashed line; ^k,2 = 0, blue-solid line for Eq. ^) from the para- 
magnetic phase (a) to the zeros of the quasi-particle energies in the 
antiferromagnetic (SDW) phase (b). Gaps open at the 'hot spots', 
that is, where the Fermi surfaces of the two flavors intersect. In this 
paper, we focus on the SDW transition that is the singular point right 
when the Fermi surfaces reconstruct. The C4 lattice symmetry of the 
original fermions is preserved. 



are 

^k,Qr ~~ ~2^Q. COS ^jk; 2.tQ; y cos ky — /Iq; . (2) 

A consistent mapping to "physical" fermions can be achieved 
with the an anisotropic choice of hoppings 1 16], t> = I and 
^< = 0.5 and: = t2,x = -t<, h^y = t2,y = -t> yielding 
the Fermi surfaces shown in Fig.[T] 

A mean-field analysis of Eq. ([T]) predicts an antiferromag- 
netic spin-density wave (SDW) ground state at r = 1.34 which 
spontaneously breaks the spin SU(2) symmetry of Eq. ([T]). 
The Fermi surface topology "reconstructs" and gaps open at 
the hot spots as shown in Fig. [T] On a mean-field level, the 
SDW transition at zero temperature of Eq. ([T]) is first order, 
as was also found in related single-band models for electronic 
antiferromagnets (171 [TSl. In the present paper, we however 
focus on continuous SDW transitions at zero temperature. 

Key results. We now describe the three key results of this 
paper, (i) We find an infrared strong-coupling fixed point for 
the Yukawa-coupling which governs the RG flow of the 
coupled Fermi-Bose action down to the lowest scales A ^ 0. 
This induces scaling relations among the anomalous expo- 
nents for the Fermi velocity, the quasi-particle weight and the 
Yukawa vertex, (ii) Both the quasi-particle weight and the 
Fermi velocity vanish as a power-law when scaling the mo- 
menta toward the hot spot; the Fermi velocity slower than the 
quasi-particle weight, (iii) The (quantum) dynamical scaling 
of the electronic single-particle and collective spin fluctua- 
tions follows from an emergent dynamical exponent attaining 
the same (fractional) value for both, fermions and boson. 

As explained in more detail in the Supplemental Material, 
our RG analysis is based on the (formally exact) flow equa- 
tion for the efl'ective action F^ ij/, ^j, the generating func- 
tional for one-particle irreducible correlation functions in the 
form derived by Wetterich |141 [ISj . The regulator R intro- 
duces a cutofl' dependence into the efl'ective action so that F^ 



smoothly interpolates between the bare action, Eq. ([T]), at the 
ultraviolet scale F^^"^"^ i//, = F^^v ^nd the fully 

renormalized efl'ective action in the limit of vanishing cutofl': 
limA^o 0] = r ij/, ^j. The Wetterich equation has 

a one-loop structure and in a vertex expansion the /^-functions 
for the n-point correlators are determined by (cutofl' deriva- 
tives of) one-particle irreducible one-loop diagrams with fully 
dressed propagators and vertices. Upon self-consistent inte- 
gration of the coupled set of /^-functions, contributions of arbi- 
trary high loop order are generated. In this paper, we truncate 
the efl'ective action to the full fermion two-point function (in- 
cluding a fermion self-energy 2^(^o. k)), the full bosonic two- 
point function (including a bosonic self-energy 2^(^o. Q)). and 
the Yukawa coupling A^. 

The centerpiece of our analysis is the flow equation for the 
Yukawa coupling: 




where (l^^ = [A^f /(A ^Jzf^Zf^ a/^/i^/2) rescaled by 

the frequency (Zfi) and momentum i^fi) derivatives of the 
fermion self-energy generated under the RG flow as per Fig. [2] 
(a). The power-law divergences as well as all other non- 
universal contributions to the flow of the two fermion self- 
energy factors and the Yukawa coupling itself are absorbed 
into the anomalous exponents: 

JlogZ^^ ^ dlogAf^ ^ dlogA^ 

"^^^^""Tl^' "^^^^ JlogA ' "^^^"""Jl^- 

(4) 

is driven by the direct contribution to the flow of 
exhibited in Fig. |2] (b). All couplings are projected to zero 
fermionic frequency, a discrete set of fermionic momenta on 
the Fermi surfaces, and zero bosonic frequency and momenta. 
This is where the most singular renormalizations occur. 




FIG. 2: Diagrammatic representation of the flow equation for the 
fermion self-energy E^(^o.k) (a) and the Yukawa coupling (b). 
Straight lines denote Fermi propagators of flavor 1 and 2, wiggly 
lines boson propagators are endowed with a regulator (specified 
in the Supplemental Material). Intersections of wiggly with straight 
lines represent the Yukawa coupling. The cutoff-derivative with re- 
spect to R^ is implicit. All propagators and vertices are "dressed" 
self-consistently and are functions of A. 
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FIG. 3: (Color online) Quantum critical RG flows of the Yukawa 
coupling and the anomalous exponents at the hot spot kns . The fixed- 
point values are = 2.38, r/z^ = 0.78, r/Af = 0.44 and //yuk = 0.11. 
The scaling plateaus for 5 > 6 depicted over ~ 4 orders of magnitude 
would be attained indefinitely but are limited by the numerics only. 
The infrared is to the right of the plot (A = Auv^~^)- 



A self-consistent numerical solution of the flow equations 
for the Yukawa vertex A^, the fermion self-energy ^^(^o^ k) 
and the boson self-energy ^^(qo, q) (to be described below) is 
attracted toward an infrared strong-coupling fixed point. As 
can be read oflT from Fig. |3j the /5-f unction for the Yukawa 
coupling, Eq. vanishes for ^ > 6 resulting in a scaling 
relation for the fermion and Yukawa anomalous exponents: 



dlogA^ 
J log A 



(8) 



where we dropped the flavor index as they become degener- 
ate at the hot spot. A similar strong-coupling fixed-point and 
scaling relations (without singular vertex corrections) have re- 
cently been obtained at the QCP of Dirac cone toy model be- 
tween a semimetal and a superfluid | 21 1. 

The numerical values of the exponents (see Fig. [3]) deter- 
mine the scaling behavior of the fermion propagator Eq. ^ 
and the associated dynamical exponent Zf. The Yukawa ver- 
tex diverges as a power-law 



Specifically, the inverse quasi-particle weight is computed 
from the flowing self-energy by 1 19] 



/A - 1 _ 

/I ~ 



d 
diko 



(5) 



where kp is a momentum on the Fermi surface and the initial 
condition is Z^^^^ = 1 . The momentum renormalization factor 
is obtained from a momentum gradient of the fermion self- 
energy 



1 + 



|nk,i-Vi:^^(/:o,k)| 



l^o=0,k=kF 



(6) 



with the initial condition A^J^^ = 1. Here, V = (dk^,dky^ and 
nk 1 is unit normal vector onto the Fermi surface of flavor 1. 
We shall see below that the momentum gradient scales dif- 
ferently than the frequency derivative at the quantum critical 
point. In a diflferent context, for Fermi systems with van Hove 
singularities, this asymmetry was established to all orders in 
perturbation theory by Feldman and Salmhofer 1201 . Neces- 
sary conditions to discover this are: (i) the co-dimension of 
the Fermi surface manifold is greater than zero (it is zero in a 
one-dimensional Fermi systems) and (ii) one includes the ad- 
ditional, relevant transversal momentum direction parallel to 
the Fermi surface into the analysis. 

With these definitions, the scale-dependent "dressed" 
fermion propagator which occurs self-consistently in all RG 
equations becomes 



= —^ = 



(7) 



with 



resembling the quasi-particle weight at 
low energies and the eflTective modulus of the Fermi velocity 



n 



1 



1 



^%uk ^0.11 



(9) 



A can be associated with the momentum distance from the hot 
spot; at A = the hot spots are resonantly connected by the 
ordering wave vector K of the incipient spin-density wave. At 
the hot spot, the fermionic quasi-particle weight vanishes as a 
power-law: 



I\^'f = A' 



0.78 



(10) 



destroying the Fermi liquid character of fermionic quasi- 
particle excitations. In a non- self consistent calculation we can 
also compute the fermion self-energy from Eq. ^ away from 
the hot spot by solving the flow equations evaluated at general 




FIG. 4: (Color online). Infrared values of the momentum resolved in- 
verse quasi-particle weights Z^:[^^[ko = 0, k^, ky] non- self consistently 
computed from Eq. ^ along the Fermi surface. Fig.[5]exhibits flows 
of the corresponding exponents for the six data points closest to the 
maximum/hot spot on the right flank. Here the hot spot is located at 
kus,y = 2.0944 and ^hs,;c = 1.0472. 
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FIG. 5: (Color online). Non-fermi liquid regimes at intermedi- 
ate scales of the anomalous exponent for the quasi-particle weight 
^Zfilko = 0,^;^,^^] for six choices of momenta progressively ap- 
proaching the hot spot (corresponding to the 6 data points closest to 
the maximum/hot spot on the right flank of Fig.|4]). The momentum 
k6 is furthest from the hot spot and ki is closest to it. The infrared is 
to the right of the plot (A = Auv^~^)- 



fermionic momenta k. The result for a momentum cut along 
the Fermi surface is exhibited in Fig.|4] The renormalization 
of the quasi-particle weight is strongly peaked around the in- 
tersection of the Fermi surfaces at the hotspot. Away from 
the hot spot, the suppression of the quasi-particle weight is 
less pronounced leading to asymptotically vanishing anoma- 
lous exponents in the infrared A ^ (Fig. [5]). Nevertheless, in 
the vicinity of the hot spot, magnetic fluctuations are still very 
strong leading to sizable non-Fermi liquid scaling regimes 
at intermediate scales with the maximum progressively ap- 
proaching the hot spot value ^/Z/J^o = = kus,x,ky = 
kns,y] = 0-78 for momenta closer to it. 

In the numerics for Fig. |4j we stopped the flow at ^ = 7 
(recall that A = Auye'^) leading to finite (but very large) val- 
ues of Zfi even at the hot spot. We used a momentum cut of 
100 points producing for each grid point in Fig. |4] the scale- 
resolved flows shown in Fig. [5] 

The Fermi velocity vanishes as well but with a smaller ex- 
ponent 



A' 



0.34 



SO that the dynamical exponent for the fermions is 
Zf = l-\-r]Zf -riAf = 1.34 . 



(11) 



(12) 



An important ingredient to the scaling laws above is the self- 
consistently flowing boson propagator 



-1 



(13) 



where is an infrared momentum regulator that vanishes 
for A ^ 0. The flow equation for the bosonic self-energy 
Sf(^o, q) is a convolution of two fermion propagators with 




FIG. 6: Particle-hole bubbles used for the flow of the boson self- 
energy in Eq. |T3] ). All propagators and vertices are "dressed" self- 
consistently and depend on A. 



two Yukawa vertices (Fig. [6]). The asymptotic static and dy- 
namic scaling of the spin fluctuation propagator is given by 

lim \D\q^, q)l"' ~ A^-^ \q^\ + ~ \q^\'-'' + q^ . (14) 

with = 0.66. Remarkably, the boson dynamical exponent 



Zb = 2-7]z, = 1.34 = z/ , 



(15) 



takes the same value as the fermion dynamical exponent. It 
is a distinguishing feature of this infrared fixed-point of elec- 
trons in metals at a spin-density wave transition that the dy- 
namical exponent attains fractional value diff'erent from 1 
(which is the exact value for quantum-critical fermion sys- 
tems with Lorentz-symmetry, see Ref. |22| and references 
therein) and diff'erent from 2 (which is the mean-field value 
of the Hertz theory iQ). Our fermion anomalous dimensions 
and z can be mapped to those of Ref. iJJj for values of the 
Fermi velocity-anisotropy in a range around a ^ 0.5, and 
upon ignoring the marginal RG flow of a (which is implicitly 
assumed in ([6])); our boson anomalous dimension renormaliz- 
ing the q^ term in the propagator is essentially zero, and we 
trace this to difl'erences in the RG scheme from Ref. \TT\ 

It would be interesting to classify all relevant operators to 
our fixed point, and investigate the stability of our strong- 
coupling fixed point further. As a first simple step in this 
direction, we have extended the truncation for the fermion 
dispersions to allow for changes in the Fermi surface curva- 
ture (keeping the position of the hot spot fixed). A scale- 
dependent that modifies the hoppings, ti^x/y ti,xiy + 5'^ 
and t2,xiy h,x/y - does the job. We found only relatively 
small, finite renormalizations of a^. However, a proper self- 
consistent investigation of a flowing Fermi surface with the 
full dispersion used in this paper requires an advanced trunca- 
tion and likely also a self-consistent determination also of the 
position of the Fermi surfaces and the hotspots as a function 
of A. Potential tendencies toward magnetic ordering at in- 
commensurate wave vectors might also be captured that way. 
Such a state-of-the-art truncation was recently presented for 
self-energy flows in the repulsive Hubbard model close to van 
Hove filling |23 |. 

Other promising future directions are the inclusion of (d- 
wave) superconductivity [24], an extension to the quantum- 
critical regime at finite temperatures, and the exploration of 
the antiferromagnetic phase with broken symmetry close to 
the quantum critical point 130 J . 



5 



Acknowledgments. We thank E. Berg, C. Honerkamp, E. 
G. Moon, and M. Punk for useful discussions. This research 
was supported by the DFG under grant Str 1176/1-1, by the 
NSF under Grant DMR- 1103860, and by the Army Research 
Office Award W91 INF- 12- 1-0227. JL is also supported by 
the STX Foundation. 



* Electronic address: |pstrack@physics.harvard.edu| 
[1] L. Taillefer, Annual Review of Condensed Matter Physics 1, 51 
(2010); N. Doiron-Leyraud and L. Taillefer, arXiv: 1204^0490 
T Helm et al , Phys. Rev. Lett. 105, 247002 (2010). 
Y. Nakai et al, Phys. Rev. Lett. 105, 107003 (2010). 
A. W. Overhauser Phys. Rev. 128, 1437 (1962). 
J. A. Hertz, Phys. Rev. B 14, 1165 (1976). 
T. Moriya, Spin Fluctuations in Itinerant Electron Magnetism, 
Springer- Verlag, Berlin (1985). 

A. J. Millis, Phys. Rev. B, 48, 7183 (1993). 
S. Sachdev, A. V. Chubukov, and A. Sokol, Phys. Rev. B 51, 
14874 (1995). 

Ar. Abanov, A. V. Chubukov, J. Schmahan, Adv. Phys. 52, 119 
(2003). 

Ar. Abanov and A. V. Chubukov, Phys. Rev. Lett. 84, 5608 
(2000); Ibid 93, 255702 (2004). 

M. A. Metlitski and S. Sachdev, Phys. Rev. B 82, 075128 
(2010). 

S.-S. Lee, Phys. Rev. B 80, 165102 (2009). 
M. A. Metlitski and S. Sachdev, Phys. Rev. B 82, 075127 
(2010). 

W. Metzner et al. Rev. Mod. Phys. 84, 299 (2012). 
J. Berges, N. Tetradis, and C. Wetterich, Physics Reports 363, 
223 (2002). 

E. Berg, M. Metlitski, and S. Sachdev,|arXiv: 1 206.0742 1(2012). 
J. Reiss, D. Rohe, and W. Metzner, Phys. Rev. B 75, 075110 

(2007) . 

B. L. Altshuler, L. B. loffe, and A. J. MilHs, Phys. Rev. B 52, 
5563 (1995). 

C. Honerkamp and M. Salmhofer, Phys. Rev. B 67, 174504 
(2003). 

J. Feldman and M. Salmhofer, Rev. Math. Phys. 20, 275 (2008). 
P Strack, S. Takei, and W. Metzner, Phys. Rev. B 81, 125103 
(2010); B. Obert, S. Takei, and W. Metzner, Ann. Phys. (Berlin) 
523, 621 (201 1). 

L. Janssen and H. Gies, arXiv: 1208.3327 (2012). 
K. U. Giering and M. Salmhofer, arXiv: 1208.6131 (2012). 
A. Sedeki, D. Bergeron, and C. Bourbonnais, Phys. Rev. B 85, 
165129 (2012). 

P Strack, R. Gersch, and W. Metzner, Phys. Rev. B 78, 014522 

(2008) . 

F. Schuetz, L. Bartosch, and P. Kopietz, Phys. Rev. B 72, 
035107 (2005). 

C. Drukier, L. Bartosch, A. Isidori, and P. Kopietz, Phys. Rev. 
B 85, 245120 (2012). 

H. Gies and C. Wetterich, Phys. Rev. D 65, 065001 (2002). 
F. Wegner, Grassmann- Variable, Lecture Notes, http : //www . 
tphys ■ uni-heidelberg . de /~wegner/VorlQ9_lQ . html 
' Universitat Heidelberg (1998). 
[30] For example, by generalizing Ref.'25'from the superfluid 0(2) 
case to the staggered 0(3) case for the spin-fermion model. 



SUPPLEMENTAL MATERIAL 

Our results are obtained from the renormalization group 
flow of the action Eq. ([T]) at the quantum-critical point (r = 0) 
under the formally exact evolution equation 1 141 



^r\x^x] = lstr[R^[Tf^\x,x]^R^]~'] • (16) 



d 

dA' 

Y^^^^ is the second derivative with respect to the fields defined 
below. is a matrix containing A-dependent cutoff' func- 
tions that regularizes the infrared singularities of the fermion 
and boson propagators. The dot is short-hand notation for a 
scale-derivative = OaR^. Both sides of this equation are 
projected onto a "super"-field basis ;f containing fermionic 
and bosonic entries: 

CPyik) 

X(k)= (17) 

and its conjugate-transposed ^(^). Str is a "super" trace over 
frequency, momenta, and internal indices and installs an addi- 
tional factor of -1 for contributions from the purely fermionic 
sector of the trace of Grassmann- valued matrices. We will 
solve Eq. ([16]) in a vertex expansion truncating any gener- 
ated vertices beyond the Yukawa vertex. The ffowing fermion 
self-energy I^^(ko, k) and the boson self-energy l^^iqo, q) are 
parametrized in a derivative expansion keeping the Fermi sur- 
faces fixed. 



The cutoff' matrix in Eq. ( 16 ) is given by 



R^ - diag(/?^'^,/?^^,7?^^,7?^[^^,7?^^^^ 

where one is in principle free to choose the fermion and boson 
cutoff' scales A^(A) and A/(A) and associated regulator func- 
tions Rhj independently 1261 [27|. The corresponding "flow 
trajectories" in cutofl' space (in the plane of Fig. [7]) from the 
bare action (red dot) to renormalized, efl'ective action (green 
dot) will be difl'erent. We will choose the trajectory along 
the arrows illustrated in Fig. |7] that is we take Af ^ 
and Rf ^ before integrating out order parameter fluctua- 
tions which are excluded for momenta smaller than A^. The 
fermions are however not discarded as in the Hertz theory 151, 
but coupled self-consistently into the flow for all A e {A^^, 0} 
thereby imposing important boundary conditions for the inte- 
gration of order parameter fluctuations down the vertical axis 
in Eig.jT] 
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FIG. 7: (Color online) Illustrative flow trajectory in cutoff space. At 
each step AA^ of the integration over bosonic momenta along the 
vertical axis the entire range of fermionic momenta is swept over 
(grey- striped box). 



For the bosons, we use a Litim cutoff for momenta 

= <; = = < = (-q2 + a2) 0(a2 - q2) , (19) 

where is bosonic momentum renormalization factor to be 
specified below. In the following, we will set A^ = A. The 



fermionic entries in Eq. ( 18 ) are zero. 

The fermionic matrix elements of the generalized matrix 

propagator [f^^^^ \x^x] + occurring in Eq. ( 16 ) become: 



fl,(T 

x=x=o 

ki=k2=k 



-1 



(20) 



and analogously for the other flavor and spin components. 
The explicitly cutoff-dependent boson propagators are 



~5 



b,x 

x=x=o 

q\=q2=q 



-1 



q2+r+2:^(^o,q) 
-1 

A2+r+2:^(^o,q) 



lql>A 
lql<A 



(21) 



and analogously for the other spin projections z. The func- 
tional derivatives are evaluated at zero fields here, as we ap- 
proach the QCP from the paramagnetic phase. 

The flow equation for the fermion self-energy (depicted di- 
agrammatically in Fig.|2]) is 

dtj!}^ [/:o, k] = 3 (^^)' r G)^{k + ^)Z)f (^) , (22) 



and similarly for flavor 2 upon interchanging 1 2. We use a 
short-hand notation encapsulating frequency, momentum in- 
tegrations and a cutoff derivative with respect to the bosonic 
cutoff function: /^,^ = ftl^ f-^^\4 

The prefactors and signs of the flow equations are computed 
by comparing coefficients between the left-hand- side and the 
right-hand- side of Eq. ([16]) as outlined in Sec. II of Ref. |28] 
The 11x11 Grassmann- valued (super-) matrices are evaluated 
using the GrassmannOps.m package in Mathematica. How to 
take a supertrace can be found in Ref. |29l 

The boson self-energy is determined self-consistently from 
the particle hole bubble (Fig. [6]) at all stages of the flow: 

S^(^o, q) = - (n^(^o, q) - n^(0, 0)) (23) 
= 2 {A'^f Jj{G%(k + q)- G%{k)) G%(k) 
^G%(k){G%(k^q)-G%(k))] 

The following ansatz captures the leading frequency and 
momentum-dependence of the particle-hole bubble: 



A^2 



mqo,q) = Z^\qo\^A^q 



(24) 



At the yellow dot in Fig. |7] the Fermi propagators are still 
Fermi-liquid like (^^^^ = 0) because we have not yet inte- 
grated out any order parameter fluctuations which, by Fig.|2j 
generate a finite fermion self-energy. At that point, the coef- 
ficients Z^^^, A^^^ take finite numerical values. At all stages 
of the flow, when integrating the flow down the vertical axis 
of Fig. |7j the bosonic Z-f actor and A-f actor are determined 
self-consistently according to the prescription: 



7A _ _ 



n^^o,o)-n^^o,o) | 

qo 1^0= A 

n^(o, q) - n^(o, 0) I 



\q^^A,qy^O 



(25) 



This allows them to pick up potentially singular renormaliza- 
tions during the flow. The boson momentum factor is isotropic 
in momentum space; interchanging q^ % delivers the same 



value for A^. 



The flow equation as per Fig. [2] for the Yukawa coupling is 

dAA^ = -U^f f G%(k^q)G%(k^q)D^,(q)\ . 

^ ^ JgR 1/^0=0, k=kHs 



(26) 



We now give the explicit expressions of the flow equations 
( |22|23|26| ). To that end, it is convenient to use the rescaled 
variables Z^ = ^k,i = ^ as well as rescaled momenta: 
^o=|,§o=f,§. = |,and§, = |. 

For the the fermionic frequency exponent, there is 
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(/§0-|l^^2l^kHs+q,2f (Z,^l§0l+A^f ' 



(27) 



and similarly (1 2) for flavor 2. The frequency integral over §o can be performed analytically so that at each step of the 
flow, two-dimensional integrations over and qy have to be performed numerically. The Yukawa anomalous exponent contains 
fermion propagators of both flavors: 



^'dq^ r^'-~^ydq, 
-1 27r j-yr^ 

For the flow of the fermionic momentum factors we use the 



1 



1 



1 



(28) 



projected kx and ky derivatives of Eq. ( 22 ) 



^% y ~ 'z^' ^i^h ^Zfi given in Eq. (|27|), these take the form 



dAA%^y = nk^,idk^dA^%[ko,k]\ 

with the initial conditions = Afi^y ■ 
face normal projector is (similarly for flavor 2) 



(29) 



li„=0,k=kHs 
<ro=0,k=kHs 

= 1. The Fermi sur- with the exponents tja 



(31) 



.11^ At 



(30) 



every step of the flow, we compute then per Eq. ([6]) 



1^/1 1 = 



The flow equations for the rescaled variables are - ^ 



|V^l,klk=kHS 

Expressions for the exponents: 



(32) 



^ > y f f An..J-i 2n J_^^ In J_ In _ |^^^^^_ (z,^l§ol + A,^)' 



^§0^ ^2f2,,sin(fcHS,y +%A) 



I 

Finally, the (rescaled) boson frequency factor and momen- turn factor are self-consistently determined from 

I 



(33) 



dk^ 
2^ 



^ /(^o + 1) - li^ji l^k,i ik^ - li^^i l^k,i 
1 



Eqs. ( |3|27|28|31|33|32|34| ) are solved numerically as a function of flow parameter 
I UV (A^^ =1). The hotspot coordinates are ^hs,x = 1.0472, ^Hs,y = 2.0944. As i 



1 



ik^ - |l3^2l^k,2 



+ (1^2) 



1 ] 


1 


ihi - ItJ^jl^k,!^ 





+ (1^2) 



(34) 



A = A^^^ ^ so that 5- = corresponds to 
initial conditions, we choose /l^^^ = 0.25, 
Z^r = = 1, and A^r = = 1. The initial values for the boson propagator are zf'^^ = 0.052 and A^"^ = 1.01 1. 



the 



/I 



/2 



